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HIGHLIGHTS 


►  Comparison  between  transport 
properties  of  Li-  and  Na-based  non- 
aqueous  electrolytes. 

►  Density,  viscosity  and  conductivity 
of  LiC104  and  NaC104  in  PC  and  yBL 
were  measured. 

►  Na-based  electrolytes  have  10-20% 
higher  conductivity  than  Li-based 
ones  at  0.5-2  M. 

►  Conductivities  were  analyzed  based 
on  a  theory  derived  from  the  pseu¬ 
dolattice  model. 
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The  solution  density,  viscosity  and  conductivity  of  MCIO4  (M  =  Li  and  Na)  solutions  in  propylene 
carbonate  and  y-butyrolactone  are  measured  at  the  concentrations  of  <1.5-2.0  mol  dnr3.  The  partial 
volume  of  the  solute,  derived  from  the  density,  of  NaC104  is  greater  than  that  of  LiC104  as  expected. 
NaC104  produces  less  viscous  and  more  conductive  solutions  than  LiC104  throughout  the  examined 
concentration  range.  Notably,  the  conductivity  of  the  NaC104  solutions  is  10—20%  higher  than  that  of  the 
LiC104  solutions  at  T/K  =  298.  The  validity  of  the  empirical  cubic  root  law,  /1(C)  =  Ap  AC1^3,  is  examined, 
where  A  and  Ap  are  the  molar  conductivities  at  the  molarity  C  and  at  infinite  dilution.  The  meaning  of  the 
slope  A  is  interpreted  in  the  theoretical  framework  of  the  pseudolattice  model. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Electrochemical  devices  based  on  nonaqueous  lithium  ion  elec¬ 
trolytes  are  now  ubiquitous  in  modern  life  from  mobile  phones  to 
electric  vehicles.  The  electrical  conductivity  and  viscosity  of  the 
electrolyte  solutions  are  a  pair  of  fundamental  transport  phenomena 
that  plays  an  important  role  in  determining  the  performance  of 
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these  electrochemical  devices.  Many  nonaqueous  solutions  doped 
with  Li  salts  have  a  maximum  specific  conductivity  <7  at  around 
1  mol  dm-3.  The  maximum  in  a  at  a  certain  concentration  can  be 
qualitatively  explained  by  the  competition  between  the  increase  in 
the  number  of  charge  carriers  and  the  dragging  (negative)  effect  of 
interionic  interactions.  While  a  number  of  experimental  data  has 
been  collected  to  cover  these  high  concentrations  [1—8],  the  theo¬ 
retical  understanding  of  the  ionic  transport  in  concentrated  solu¬ 
tions  is  not  yet  fully  established  (see  Section  3.3). 

Meanwhile,  attempts  to  replace  Li  with  Na  as  a  charge  carrier  in 
the  electrochemical  devices  have  been  recently  drawing  attention 
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Table  1 

Symbols  in  the  present  study. 

Symbol  Quantity 

C  Molarity  (mol  concentration) 

C*  Molarity  at  which  a  is  maximum 

i)  Viscosity  of  solution 

T]0  Viscosity  of  pure  solvent 

yr  Relative  viscosity,  17/ijo 

A  Molar  conductivity.  <r/C 

A’  OnnJC 

M]  Molecular  weight  of  solvent 

M2  Molecular  weight  of  solute 

nj  Molar  number  of  solvent 

n2  Molar  number  of  solute 

a  Specific  conductivity 

(tmax  Maximum  in  a 

p  Density  of  solution 

Po  Density  of  pure  solvent 

yd  Molar  volume  of  pure  solvent 

yd  Molar  volume  of  pure  solute  (crystal) 

yj  Partial  molar  volume  of  solute  in  solution 


due  to  the  possible  scarcity  and/or  the  geopolitical  volatility  of  Li 
resources  [9-17],  Based  on  the  advent  of  Na  ion  devices  in  the 
future,  one  finds  that  fundamental  data  on  Na-based  electrolytes 
are  lacking  in  the  literature  compared  to  the  vast  compilation  of 
information  about  the  Li-based  ones. 

In  the  present  study,  we  focused  our  attention  on  the  differences 
in  the  transport  properties  between  Li-  and  Na-salts  by  choosing 
LiClC>4  and  NaCICU  dissolved  in  two  nonaqueous  solvents, 
propylene  carbonate  (PC)  and  y-butyrolactone  (yBL)  as  model 
systems.  The  conductivities  and  viscosities  of  the  LiClC>4  and  NaC104 
solutions  are  compared  on  the  basis  of  empirical  formulas.  In 
addition,  the  conductivity  is  analyzed  in  the  framework  of  a  theory 
derived  from  the  pseudolattice  model. 

The  symbols  and  basic  properties  of  the  solvents  and  solutes 
used  in  the  present  study  are  summarized  in  Tables  1  and  2. 

2.  Experimental 


1.7  mol  dm-3  of  NaClC>4-PC,  in  which  the  concentration  of  the  latter 
one  is  limited  by  the  solubility.  The  stock  solutions  were  dried  by 
using  the  type-5A  molecular  sieves  for  72  h.  The  concentration  of 
each  electrolyte  solution  for  the  measurement  was  adjusted  by 
diluting  the  stock  solution  with  pure  solvent.  The  water  content  of 
the  solutions  was  maintained  at  less  than  50  ppm,  which  was 
determined  by  the  Karl  Fischer  titration  method  (831  KF  Coulo- 
meter,  Metrohm). 

The  temperature  of  the  solutions  was  controlled  in  a  thermostat 
within  A T/K  =  ±0.3.  The  temperature  of  the  thermostat  was  cali¬ 
brated  by  a  certified  standard  thermometer.  The  solutions  were  left 
in  the  thermostat  for  30  min  before  the  measurement  to  equilibrate 
the  temperature.  For  all  the  measurements,  three  batches  of  solu¬ 
tions  were  prepared  at  each  concentration.  The  results  were  pre¬ 
sented  as  the  mean  of  three  runs  using  each  solution  batch. 

The  density  of  the  solutions  was  measured  by  pycnometers  with 
the  nominal  internal  volume  of  10  cm3  of  which  the  accurate  value 
was  calibrated  by  measuring  the  densities  of  degassed  distilled 
water.  The  mass  of  the  fluid  was  determined  by  subtracting  the 
mass  of  the  dried  empty  pycnometer  from  that  of  the  filled  one. 

The  viscosity  was  measured  at  T/K  =  278,  288  and  298  by  the 
rotating  spindle  method  (Brookfield  DV-II±)  equipped  with  a  UL 
attachment  to  cover  the  viscosity  range  of  the  present  study.  The 
viscometer  was  calibrated  using  standard  reference  oils  (Nippon 
Grease,  JS2.5,  JS5  and  JS50)  covering  the  viscosity  range  from  1.9  to 
42  mPa  s  at  T/K  =  298. 

The  conductivities  of  the  solutions  were  measured  at  T/I<  =  288, 
298  and  308.  A  beaker-type  cell  with  four  platinum  electrodes  was 
built  in  house  [18].  The  AC  four-terminal  method  was  adopted  using 
an  electrochemical  analyzer  (IviumStat).  The  applied  current 
amplitude  was  0.1  mA  under  which  the  amplitude  of  the  voltage 
response  ranged  from  ca.  1  to  1000  mV  depending  on  the  solution. 
The  resistance  was  calculated  from  the  I—V  relation  at  20  kFIz  at 
which  the  observed  conductivity  was  confirmed  to  be  independent 
of  the  frequency  and  free  of  noise.  The  resistance  was  converted  to 
the  specific  conductivity  a  by  calibrating  the  cell  with  0.01  and 
0.1  mol  dm  3  KC1  aqueous  solutions  of  which  the  specific  conduc¬ 
tivities  are  known  to  be  1.4085  and  12.853  mS  cm1,  respectively,  at 
T/K  =298  [19], 


Preparation  of  the  solutions  and  the  setup  of  the  cells  for  all  the 
measurement  were  carried  out  in  a  dry-box  in  which  the  dew  point 
of  air  was  maintained  below  210  K.  Lithium  perchlorate  (IJCIO4, 
nominal  purity  97%)  and  sodium  perchlorate  (NaCKXi,  96%)  were 
purchased  from  Kanto  Chemical.  The  solutes  were  used  without 
further  purification,  as  we  confirmed  that  the  solute  recrystallized 
from  tetrahydrofuran  showed  no  significant  difference  in  conduc¬ 
tivity  from  the  one  without  recrystallization.  The  purities  were 
taken  into  account  in  weighing  the  solute.  Solvents  were  purchased 
from  Kishida  Chemical;  propylene  carbonate  (PC,  purity  >  99.5%, 
water  content  <  30  ppm)  and  y-butyrolactone  (yBL,  purity  >  99.5%, 
water  content  <  30  ppm).  We  first  prepared  the  stock  solutions  of 
2.0  mol  dm-3  of  LiC104-PC,  LiC104-yBL  and  NaC104-yBL;  and 


Basic  properties  of  the  solvents  and  solutes. 


Solvent 

PC  CH3(C2H3)02C=0 

yBL  (CH2)3OC=0 

Molecular  weight,  M,/ g  mol"1 

102.1 

86.1 

Relative  permittivity,  er 

64.98 

41.78 

Dipole  moment 

4.94D 

4.12D 

Molar  volume,  vf/cm3  mol-1 

85.1 

76.6 

Solute 

LiC104 

NaC104 

Molecular  weight,  M2/g  mol"1 

106.4 

122.4 

Molar  volume3,  V|/cm3  mol-1 

43.8 

49.2 

d  In  crystalline  state  calculated  from  the  crystallographic  data. 


3.  Theoretical  frameworks 

3.1.  Solution  density  and  partial  volume  of  the  solute 


The  partial  volume  of  the  solute  V2,  i.e.,  the  increase  in  the 
solution  volume  when  an  infinitesimal  amount  of  the  solute  is 
added  to  the  solution  at  a  given  concentration,  is  defined  by 


where  V  is  the  volume  of  the  solution,  and  ni  and  n2  are  the  molar 
numbers  of  the  solvent  and  solute,  respectively.  Note  that  V2 
includes  not  only  the  volume  of  the  solute  molecule  but  also  the 
change  in  the  volume  occupied  by  the  solvent  molecules  that 
solvate  the  solute  molecule.  If  the  solution  density  is  known  as 
a  function  of  the  concentration  p  =  p(C),  one  can  derive  the  partial 
molar  volume  as  follows.  The  volume  of  the  solution  is  related  to 
the  density  as 

v  =  n-iMi  +  n2M2 


where  Mi  and  M2  are  the  molecular  weights  of  the  solvent  and 
solute.  For  all  the  systems  examined  in  the  present  study,  as  shown 
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later  in  Section  4.1,  the  density  p  linearly  varies  with  the  molarity  C 


P  = 


Po  +  aC 


(3) 


(4) 


where  po  is  the  density  of  the  pure  solvent  and  a  is  a  constant. 
Introducing  Eq.  (4)  into  Eq.  (2)  and  solving  for  V ,  one  gets 


v  =  niM,  +n2(M2-a) 
Po 


=  mv?  +  n2 


M2  -  a 
Po 


(5) 


where  1/°  is  the  molar  volume  of  the  pure  solvent.  Hence,  the 
partial  molar  volume  of  the  solute  Eq.  (1)  in  this  case  is  given  by 


V2 


M2-  a 
Po 


(6) 


Namely,  the  linear  dependency  of  the  solution  density  on  the 
molarity,  Eq.  (3),  suggests  the  constant  partial  volume  of  the  solute 
independent  of  the  concentration. 


3.2.  Viscosity 


The  Jones— Dole  empirical  relation  between  the  viscosity  p  and 
the  molarity  C  is 

Vr=—  =  1  +AC1/2  +BC  +  DC2  (7) 

Vo 


where  p  and  po  are  the  viscosities  of  the  solution  and  the  pure 
solvent,  respectively,  and  A,  B  and  D  are  parameters  [20],  The /l-  and 
B-terms  can  be  physically  interpretable  [21  —25].  On  the  other  hand, 
although  the  interionic  interactions  seem  to  be  related  to  it  [26],  the 
parameter  D  is  currently  considered  empirical.  The  71-term  is 
significant  only  at  low  concentrations,  and  thus  often  omitted  at 
higher  concentrations,  e.g.,  at  C/ mol  dm-3  >  0.05.  In  such  case. 


which  the  well-known  Debye— Huckel— Onsager  theory  provides 
the  theoretical  basis  (See  Section  3.3.3).  However,  at  higher 
concentrations,  e.g.,  0.1  <  C/mol  dm-3  <  2  at  which  many  electro¬ 
chemical  devices  work,  Eq.  (10)  with  n  =  3  is  empirically  known  to 
better  represent  the  experimental  data.  If  n  =  3  in  Eq.  (10), 


1  =  Aq-AC1^, 

and  the  specific  conductivity  a  is  given  by 

(11) 

r=AC  =  C(A0-AC1/3'). 

(12) 

Solving  dff/dC  =  0,  one  finds  a  maximum  in  Eq.  (12)  at 


The  maximum  in  the  specific  conductivity  <rmax  at  this  molarity 
C*  is  thus 


ffmax  =  C*(A0-AC*V3) 
—  r*Ao 


(14) 


Eq.  (14)  suggests  that,  if  the  cubic  root  law  Eq.  (11)  holds,  the 
molar  conductivity  A*  at  C*,  at  which  the  specific  conductivity 
shows  the  maximum  omax,  is  one  quarter  of  the  molar 
conductivity  at  the  infinite  dilution  regardless  of  the  system; 
that  is, 


When  we  normalize  o  and  C  by  amax  and  C*  as 
o 

Cmax 


(15) 


(16) 


(17) 


Vr 


1  +BC+  DC" 


(8)  Eq.  (12)  is  reduced  to 


where  the  square  term  in  Eq.  (7)  is  replaced  by  an  n-th  power  term 
to  deal  with  the  general  case  in  the  discussion  below.  Hence, 

=  B  +  DCn  \  (9) 

Namely,  the  plot  (ijr  —  1  )/C  versus  Cn_1  should  give  a  straight  line 
of  which  the  y-intercept  and  the  slope  give  B  and  D,  respectively. 
We  will  discuss  the  validities  of  these  empirical  equations  in 
Section  4.2. 


y  =  4x(l-|x^.  (18) 

Notice  that  Eq.  (18)  no  longer  includes  any  adjustable  parame¬ 
ters  that  depend  on  the  system;  i.e.,  the  law  of  corresponding  state. 

3.3.2.  Walden’s  rule 

Based  on  the  continuum  model,  a  rough  relation  between  the 
molar  conductivity  A  and  the  viscosity  p  is  given  by  the  Walden  rule 
which  is,  for  a  singly  charged  species,  derived  from 


3.3.  Conductivity 

Classical  theories  of  electrolyte  conductivity  are  well  docu¬ 
mented  in  many  textbooks  such  as  Refs.  [27,28].  In  this  section,  only 
an  outline  that  is  pertinent  to  the  present  discussion  is  described. 

3.3.1.  Empirical  cubic  root  (C^3)  law  of  the  molar  conductivity 
Suppose  the  molar  conductivity  A  obeys  the  following  formula 

A  =  A0  -  AC1!"  (10) 


(19) 


where  kB  is  the  Boltzmann  constant;  F  the  Faraday  constant;  e  the 
unit  charge;  and  D,  r  and  X  are,  respectively,  the  diffusion  coeffi¬ 
cient,  the  radius  and  the  ionic  conductivity  of  the  drifting  entity. 
This  relation  gives  the  Walden  product  as 


(20) 


where  Ao  is  the  molar  conductivity  at  infinite  dilution.  In  the  vicinity 
of  infinite  dilution,  Kohlrausch  experimentally  observed  n  =  2  to 


(21) 
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or  in  terms  of  the  molar  conductivity  A, 


^ocl  +  1  (22) 

where  r+  and  r  are  the  radii  of  the  cationic  and  anionic  species. 
Because  lime  >o?)  =  770,  i.e.,  independent  of  the  solute,  one  gets 

4^+^  (23) 


for  a  given  solvent  system. 

3.3.3.  Pseudolattice  model 

Since  the  advent  of  the  Debye-Hiickel-Onsager  (DHO)  theory, 
the  model  has  undergone  many  modifications  to  extend  its  appli¬ 
cable  range  of  concentrations.  Several  other  approaches  outside  of 
the  DHO  theory  have  also  been  formulated,  such  as  the  mean 
spherical  approximation,  hyper-netted  chain  model,  etc.  (See,  for 
example  Ref.  [28]).  These  theories,  however,  have  not  yet  succeeded 
in  satisfactorily  explaining  the  general  behavior  of  the  conductivity 
at  the  molarity  around  C/ mol  dm3  ~  1.0  or  higher.  The  Casteel- 
Amis  formula,  an  entirely  empirical  equation  of  which  the  fitting 
parameters  have  no  physical  meaning,  is  widely  used  to  precisely 
represent  the  experimental  data  at  high  concentrations  [29], 
Recently,  a  phenomenological  treatment  in  the  framework  of  the 
dissipation-feedback  scheme  [30],  or  theories  based  on  the  pseu¬ 
dolattice  model  [31—33]  have  been  successfully  applied  to  the 
conductivity  of  the  highly  concentrated  electrolytes  or,  in  some 
cases,  of  ionic  liquids.  (See  these  references  and  those  cited  therein 
for  more  details  as  well  as  the  recent  development  of  other 
theories.) 

In  the  present  study,  for  the  sake  of  a  straightforward  compar¬ 
ison  between  the  Li-  and  Na-salts,  we  adopted  the  formalism 
applied  by  Chagnes  et  al.  [31,32],  The  pseudolattice  model  has 
emerged  through  the  observation  that  certain  properties  of 
concentrated  ionic  solutions  can  be  well  explained  by  assuming 
that  the  solution  is  composed  of  a  hypothetical  lattice-like  structure 
[34,35],  A  modification  of  the  DHO  theory  by  the  pseudolattice 
model  was  proposed  as  follows:  In  the  DHO  theory,  the  molar 
conductivity  A  is  expressed  as 

A  =  Ag  —  {A\  +  A2Ao)k  (24) 

where  k-1  is  the  mean  radius  of  the  ionic  cloud  around  the  central 
ion  under  consideration.  Ai  and  A2  are  respectively  related  to  the 
electrophoresis  and  relaxation  of  the  ionic  cloud  through  the  drift 
of  the  central  ion.  For  the  l:l-electrolyte, 


Ai 


NAe2 

3itt] 


(25) 


_  /N  e2  V 
“  \y  ^rNAkBTj 


However,  at  high  concentrations  like  C/ mol  dm3  ~  1,  the  radius 
of  the  ionic  cloud  /c_1/nm  ~  0.3  which  is  of  the  same  order  of 
magnitude  as  the  size  of  the  ion  itself,  suggesting  that  the  premises 
lying  behind  the  DHO  theory  are  no  longer  valid,  nor  are  the  terms 
derived  from  the  theory,  Ai  or  A2. 

The  idea  of  salvaging  the  DHO  theory  by  the  pseudolattice 
model  at  high  concentrations  is  to  replace  the  reciprocal  radius  of 
the  ionic  cloud  k,  Eq.  (27),  by  the  following  one  [38].  The  lattice 
energy  per  ion,  £L,  for  a  l:l-valent  ionic  crystal  is 


where  a  is  the  Madelung  constant  and  ro  is  the  nearest  neighbor 
distance  in  the  lattice.  In  the  pseudolattice  model  of  the  electrolyte 
solution,  the  ions  in  the  electrolyte  solution  are  assumed  to 
constitute  a  loose  lattice-like  structure  of  which  the  nearest 
neighbor  distance  is  equal  to  the  mean  distance  between  cation  and 
anion;  i.e., 


where  V  is  the  volume  of  the  system  and  N  is  the  total  number  of 
cations  and  anions.  The  model  identifies  the  radius  of  the  ionic 
cloud  K-1  to  be  the  one  that  gives  the  same  lattice  energy,  Eq.  (28), 
when  this  countercharged  ionic  cloud  interacts  with  the  central 


£l=pt.  (30) 

Eqs.  (28)— (30)  yield 


to 

(3,) 

=  a(2  NaQ? 

where  the  factor  “2”  appears  because  n  moles  of  a  l:l-electrolyte 
contain  a  total  of  2 n  moles  of  ions.  Notice  that  if  C  is  expressed  in 
the  unit  of  “mol  dm-3”  instead  of  “mol  m-3”,  C  in  Eq.  (31 )  must  be 
multiplied  by  103  to  conform  to  the  MKSA  system. 

If  we  substitute  Eq.  (31)  for  k  in  Eq.  (24),  the  expression  of  the 
molar  conductivity  A  based  on  the  pseudolattice  model  is  given  by 


A2 


1  e2 

2+VS127rr0ErkB7' 


(26) 


A  =  Ag  —  a(Bi  +  B2A0)O 
where 


(32) 


where  JVA  is  Avogadro’s  number;  r0  and  tr  are  the  permittivity  of  23  4  1 

a  vacuum  and  the  relative  permittivity.1  The  DHO  theory  suggests  Bi  =3^wfe2  -  (33) 


1  The  denominator  in  A2  may  look  different  by  the  factor  4itf0  from  some 
references  because  of  the  conversion  of  the  dated  electromagnetic  unit  system  into 
the  current  MKSA  system.  For  instance,  D  in  Refs.  [36,37]  and  c  in  Refs.  [27,28] 
should  read  47t£0£r  in  the  present  case. 


b2 


23  N\e2  1 

12rc(2+yf\  £ofeB  trT 


(34) 
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Molarity,  C/mol-dm'3 


Fig.  X.  Solution  density  p  at  T/K  =  298. 


Vj(in  PC)  -  V2(in  yBL)~2  -  3  cm3  mor1 
for  each  solute. 

The  larger  ionic  radius  of  Na+  than  Li+,  and  accordingly  the 
larger  size  in  the  “static”  solvation  shell  of  the  Na-salt  must  be  re¬ 
flected  in  (i).  Item  (ii)  suggests  that  yBL  more  densely  solvates  the 
solute  than  PC,  whichever  solute  is  solvated  by  them.  The  simpler 
molecular  structure  of  yBL  than  PC  may  ease  conforming  to  the 
solvation  shell. 

The  ionic  conductivity  of  Na+  at  infinite  dilution,  Aq  (Na+),  is 
higher  than  that  of  Li+,  Ao  (Li+)  [39].  As  Eq.  (23)  suggests,  the  radius 
of  the  drifting  cationic  agency,  r+,  is  considered  larger  for  Li+  than 
for  Na+  because  the  former  more  strongly  attracts  the  solvent  than 
the  latter.  The  present  result  confirms  that  the  order  of  the  static 
volume  observed  in  V2  is  contrary  to  the  dynamic  volume  reflected 
in  1/A0. 


Introducing  the  numerical  values  of  the  physical  constants  and 
taking  care  of  the  compatibility  of  the  units,  one  gets  the  expression 
with  the  explicit  units  as 


A  =  Ay 

S  cm2  mol  1  S  cm2  mol 


#17451  1735  8  A0  \  f  C 

VTj/mPa  s  +  £r (T/K)  $  cm2  mo[- V  Vmol  drrT3/ 


(35) 


4.  Result  and  discussion 

Numerical  data  at  T/K  =  298  presented  in  the  present  study  as 
well  as  those  at  the  other  temperatures  are  available  in  the 
Supplementary  material. 

4.1.  Density  and  partial  volume  of  the  solute 

Fig.  1  plots  the  molarity  dependence  of  the  solution  density 
p  =  p(C)  at  T/K  =  298.  Throughout  the  measured  concentration 
range,  the  densities  are  well  represented  by  straight  lines  as  in  Eq. 
(3),  or  by  explicitly  introducing  the  units, 

=  _ft_  +  10-3 - “ - (36) 

g cm  J  gem  J  g  mol  1  mol  dm  3 

The  result  of  the  numerical  optimization  to  Eq.  (36)  is  summa¬ 
rized  in  Table  3.  Also  listed  in  the  table  is  the  partial  molar  volume 
of  the  solute  Vf  calculated  by  Eq.  (6).  One  can  find  two  features  in 
the  partial  volume: 

(i)  NaC104  occupies  ca.  3  cm3  mol-1  more  volume  than  I2CIO4 
regardless  of  the  solvent,  i.e., 

V^(NaC104)  -  V2(LiC104)~3  cm3  mol1 
for  each  solvent. 

(ii)  Either  solute  occupies  2-3  cm3  mol  'of  more  volume  in  PC 
than  in  yBL,  i.e.. 


Table  3 

Linear  relation  Eq.  (36)  between  the  density  p  and  the  molarity  C,  and  the  partial 
molar  volume  of  the  solute  V2  at  T/K  =  298. 


PC 

yBL 

UCIO4 

NaC104 

UCIO4 

NaC104 

Po/gcm  3 

1.198 

1.125 

a/g  mol1 

58.0 

70.6 

64.1 

76.2 

V2/cm3mor1 

40.4 

43.2 

37.6 

41.1 

4.2.  Viscosity 

The  viscosity  p  as  a  function  of  the  molarity  C  at  T/K  =  298  is 
shown  in  Fig.  2(a).  Either  in  PC  or  in  yBL,  LiC104  produces  a  more 
viscous  solution  than  NaC104  at  a  given  concentration.  The 
difference  between  the  LiC104-  and  NaCKXi-solutions  expands  at 
higher  concentrations.  In  yBL,  for  example,  while  the  LiCICU 
solution  is  about  10%  more  viscous  than  the  NaC104  solution  at 
C/ mol  dm3  =  1.0,  the  difference  is  more  than  30%  at 
C/mol  dm  3  =  2.0.  In  PC,  the  case  is  more  prominent  in  that  the  gap 
is  20%  at  C/mol  dm-3  =  1.0,  which  increases  to  45%  at  C/ 
mol  dm-3  =  1.5. 


Fig.  2.  (a)  Viscosity  tj  as  a  function  of  the  molarity  C  at  T/K  =  298.  (b)  (%-l)/C  vs.  C2 
(the  inset  is  vs.  C).  Solid  lines  indicate  the  function  optimized  to  Eq.  (37). 
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Table  4 

Viscosity  of  pure  solvent  tj0,  and  the  optimized  parameters  B'  and  D'  in  Eq.  (37)  at 
T/K  =  298. 

t]o  mPa  s  LiC104  NaC104 

B'  dm3  mol-1  D'  dm9  mol-3  B'  dm3  mol-1  D'  dm9  mol-3 
PC  233  1.077  1.293  1.123  0.636 

yBL  1.76  0.880  0.638  0.861  0.403 


We  found  that  the  original  Jones— Dole  equation,  i.e.,  n  =  2  in  Eq. 
(8),  does  not  well  represent  the  experimental  data  at  higher 
concentrations,  such  as  at  C/ mol  dm'3  >  1.5.  By  replacing  the  square 
term  n  =  2  with  a  cubic  term  n  -  3,  we  obtained  a  better  represen¬ 
tation  of  the  experimental  data  without  increasing  the  number  of 
fitting  parameters  as  indicated  by  the  solid  lines  in  Fig.  2(a)  and  (b): 

Vr  =  —  =  l+B'C  +  D'C3  (37) 

Vo 

or 

B'  +  D'C2  (38) 

where  we  replaced  the  symbols  B  and  D  with  B'  and  D'  to  distin¬ 
guish  from  the  original  Jones— Dole  equation  Eq.  (7).  As  shown  in 
Fig.  2(b),  the  plot  of  (?jr  -  1  )/C  versus  C2  is  closer  to  linear  than  that 
of  versus  C,  which  is  obviously  concave  upward  as  in  the  inset 
figure.  Table  4  summarizes  the  optimized  parameters  of  B ’  and  D'  as 
well  as  the  experimentally  measured  viscosities  of  the  pure 
solvents  i/o-  Roughly  speaking,  in  a  given  solvent,  whereas  the 
B'-parameters  of  UCIO4  and  NaC104  are  close,  D'  of  LiC104  is  greater 
than  that  of  NaC104,  which  is  responsible  for  the  increasing  gap  in 
the  viscosity  between  these  salts  at  higher  concentrations.  In  the 
present  study,  we  present  these  parameters  purely  from  an 
empirical  basis.  Although  the  original  D-parameter  is  considered 
empirical  in  the  first  place,  the  physical  interpretation  of  the 
B-term,  see  for  example  Ref  [25],  does  not  hold  any  more  for  B'. 

4.3.  Conductivity 

Fig.  3  plots  the  specific  conductivity  o  versus  the  molarity  C  at  T/ 
K  —  298.  Scales  above  the  figure  indicate  the  solvent  to  solute  molar 
ratio  ni/n2  for  each  solvent  system.2  The  data  are  well  represented 
by  the  empirical  Casteel— Amis  equation  [29]  as  indicated  by  the 
solid  line  (See  Supplementary  material).  The  specific  conductivity  a 
of  the  LiC104  solution  in  PC  at  1.0  mol  dm'3  and  298  K  is  in  our 
study  5.01  mS  cm-1,  which  falls  within  the  rather  scattered  values 
found  in  preceding  studies;  i.e.,  5.6  [1],  3.9  [41]  and  4.356  mS  cm'1 
[42],  Our  result  of  <r(LiC104in  yBL,  1.0  mol  dm'3,  298  K )/ 
mS  cm'1  =  10.4  well  coincides  with  the  one  determined  by  Chagnes 
et  al.,  10.9  mS  cm'1  [31  ].  For  the  NaC104  solution  in  PC  at  298  K,  our 
result  of  5.68  mS  cm'1  at  0.9  mol  dm'3  (=0.782  mol  kg'1)  is 
comparable  to  6.13  mS  cm'1  which  is  reported  as  the  conductivity 
of  the  same  system  at  0.78964  mol  kg'1  [3], 

The  yBL  solutions  have  higher  conductivities  than  the  PC  solu¬ 
tions,  for  which  the  lower  viscosity  of  yBL  than  PC  is  considered  to 
some  extent  responsible.  The  conductivities  of  the  NaC104  solutions 
are  always  higher  than  those  of  the  corresponding  LiClC>4  solutions. 
Especially,  in  the  concentration  range  in  which  many  electro¬ 
chemical  devices  function  (0.5  <  C/mol  dm'3  <  1.5),  the 


2  Exactly  speaking,  ti]/n2  depends  on  not  only  the  solvent,  but  also  the  solute. 
However,  the  difference  in  ni/n2  between  LiC104  and  NaC104  is  not  recognizable  in 


Fig.  3.  Specific  conductivity  1 7  plotted  vs.  the  molarity  C  of  LiC104  and  NaC104  dissolved 
in  PC  and  yBL  at  T/K  =  298.  Solid  lines  are  the  result  of  the  numerical  optimization  to 
the  empirical  Casteel-Amis  function  [29],  Scales  above  the  figure  indicate  the  solvent 
to  solute  molar  ratio  nt/n2  for  each  solvent  system. 

conductivity  of  the  NaC104  solution  is  10-20%  greater  than  that  of 
the  LiC104  solution.  So  long  as  the  conductivity  of  the  electrolyte  is 
concerned,  the  sodium-based  electrolyte  solution  can  rival  or 
surpass  the  lithium-based  counterpart. 

In  yBL,  the  Walden  product  At),  Eq.  (22),  of  the  NaC104  solution  is 
still  greater  than  that  of  the  LiC104  solution  even  at  high  concen¬ 
trations,  similar  to  the  case  of  the  dilute  solutions.  In  PC,  on  the 
other  hand,  At]  of  the  NaCICU  solution  is  not  always  higher  than  that 
of  the  LiClC>4  solution  depending  on  the  concentration.  This 
suggests  that  the  continuum  model  based  on  which  the  Walden 
rule  is  derived  is  no  longer  valid  at  high  concentrations  as  pointed 
out  by  Brouillette  et  al.  [8], 

Table  5  summarizes  the  maximum  in  the  specific  conductivity 
<7max  at  the  molarity  C*.  and  the  molar  conductivity  A*  at  C*.  i.e., 
A*  =  <Tmax/C*.  Also  listed  in  the  table  is  the  molar  conductivity  at 
infinite  dilution  divided  by  four,  Aq/4,  calculated  from  the  data 
given  in  Refs.  [39,40],  Eq.  (15)  holds  fairly  well  for  all  the  investi¬ 
gated  systems,  suggesting  the  empirical  cubic  root  law,  Eq.  (11), 
well  represents  the  solution  conductivity  at  least  in  the  vicinity  of 

The  validity  of  the  cubic  root  law  is  also  examined  in  Fig.  4  in 
which  the  molar  conductivity  A  is  plotted  versus  C1^3.  In  yBL,  the 
relation  A  =  /1(CI/3)  falls  on  a  straight  line  in  the  range  0.45  <  C1/3/ 
mol 1,3  dm'1  <  1.1  or  0.1  <  C/mol  dm'3  <  1.3.  On  the  other  hand,  the 
linearity  is  fair,  at  best,  for  the  PC  systems.  (This  issue  is  revisited 
below  in  the  discussion  on  the  law  of  the  corresponding  state.)  The 
results  of  the  numerical  optimization  of  the  experimentally  ob¬ 
tained  A  to  Eq.  (11)  are  listed  in  Table  6.  The  calculated  molar 
conductivities  at  infinite  dilution,  A0,  well  coincide  with  those  re¬ 
ported  in  Ref.  [39]  for  LiC104  and  NaC104,  and  in  Ref.  [40]  for  LiC104. 
Either  in  PC  or  in  yBL,  the  slope  of  the  straight  line  in  A  vs.  C1^3  is 
more  negative  for  NaC104  than  for  UQO4;  i.e., 

A(NaC104)  >  A(LiC104).  (39) 

Nonetheless,  the  conductivities  of  the  NaCICU  solutions  exceed 
those  of  the  LiC104  solutions  because  A(NaC104)  is  not  very  high  to 
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Table  5 

Maximum  in  the  specific  conductivity  <rmax  at  the  molarity  C*  at  T/K  =  298.  The  molar 
conductivity  A *  at  this  molarity  C*  is  compared  with  the  molar  conductivity  at 
infinite  dilution  divided  by  four,  A01 4  (cf.  Eq.  (15)). 


PC  yBL 


UCIO4  NaC104  UCIO4  NaC104 

C*/mol  dm-3 

0.8 

0.8 

1.0 

i.i 

ffmax/mS  cm-1 

5.23 

5.73 

10.4 

12.3 

a*/s  cm2  mor1 

6.54 

7.16 

10.4 

11.2 

[AotAflS  cm2  mor1 

6.69 

6.97 

9.80 

10.0 

Mo/4)b/S  cm2  mor1 

6.84 

- 

10.6 

- 

a  Calculated  from  /10  in 
b  Calculated  from  /10  in 

Refs.  [39], 
Refs.  [40], 

cancel  out,  at  higher  concentrations,  the  higher  molar  conductivity 
at  infinite  dilution,  ^(NaClCU)  >  /1o(LiC104). 

The  difference  in  the  slope  A  can  be  semi-quantitatively 
explained  by  the  pseudolattice  model.  As  Eq.  (35)  suggests,  the 
higher  Aq  leads  to  the  higher  A,  if  t]  and  er  are  assumed  those  of  the 
solvent,  not  of  the  solution,  though  this  assumption  is  arguable.  In 
effect,  Eq.  (35)  yields  the  slope  Aca i  in  Table  6,  if  we  hypothesize  the 
Madelung  constant  of  the  pseudolattice  to  be  of  the  regular  fee 
lattice,  a  ss  1.748  and  Aq  to  be  the  one  that  is  obtained  by  the  linear 
optimization  in  the  present  study.  Regardless  of  the  solute,  NaC104 
or  LiC104,  the  calculated  slope,  Acai,  explains  ca.  70%  and  90%  of  the 
slope  of  the  experimental  data  Aexp  in  PC  and  yBL,  respectively, 
suggesting  that  the  tendency,  Eq.  (39),  is  consistent  with  the  model. 
To  change  the  viewpoint,  acai  in  Table  6  is  obtained  if  we  assume, 
conversely  to  the  fee  lattice  assumption,  the  Madelung  constant 
a  to  be  a  fitting  parameter  to  reproduce  the  experimental  data. 
Whereas  aCai  strongly  depends  on  the  solvent,  it  is  not  very 


Result  of  the  numerical  optimization  of  a  (at  C/mol  dm  3  >  0.1 )  to  the  cubic  root  law 
Eq.  (11)  at  T/K  =  298,  compared  with  the  parameters  calculated  based  on  the 
pseudolattice  model. 


PC 

yBL 

UC104 

NaC104 

UCIO4 

NaC104 

AolS  cm2  mol-1 

28.0 

29.9 

41.8 

44.5 

Aexp/102  S  cm3  r 

nor4/3  2.25 

2.39 

3.11 

3.21 

Acai/102  S  cm3  n 

lol-4'3  1.64 

1.67 

2.75 

2.82 

Acal/Aexp 

0.731 

0.702 

0.886 

0.879 

«cal 

2.39 

2.49 

1.97 

1.99 

40,Aexp:  Optimize 

d  parameters  in  Eq.  ( 1 1 ) 

using  experimental  data;  Aal 

:  Calculated 

slope  in  Eq.  (35)  u 

nder  the  assumption  a  = 

=  1.748  (fee);  ( 

Xcai:  Calculated  t 

1  in  Eq.  (35) 

in  order  to  reproduce  Ae 


Fig.  5.  Normalized  conductivity  y  =  <t/em ^  vs.  normalized  molarity  x  =  C/C.  Chagnes 
et  al.:  Refs.  [31,32]  (See  also  Acknowledgement).  Solid  line  indicates  Eq.  (18). 


dependent  on  the  solute.  Using  Raman  spectroscopy,  Chagnes  et  al. 
discussed  that  the  higher  slope  of  Aexp  than  Acai  for  LiCKXpn  yBL 
might  be  related  to  their  observation  that  the  degree  of  ion  pairing, 
which  reduces  the  fraction  of  charged  species  that  contribute  to  the 
conductivity,  was  roughly  constant  over  a  wide  concentration  range 
[32],  Similarly,  the  formation  of  the  ion  pairs  and  the  solvent- 
shared  ion  pairs,  Li+  - solvent- -  C104-,  is  suggested  in  the  concen¬ 
trated  solution  of  LiC104  [43,44],  The  present  study  implies  that  the 
influence  of  the  ion  pairing  on  the  slope,  A,  is  more  prominent  in  PC 
than  in  yBL,  but  the  difference  between  LiC104  and  NaC104  is  less 
significant  than  the  difference  in  the  solvent. 

Returning  to  the  empirical  cubic  root  law,  Fig.  5  plots  the 
normalized  specific  conductivity  y  =  o7<rmax  versus  the  normalized 
molarity  x  =  C/C*,  in  which  (jmax  and  C*  were  calculated  from 
A0  and  Aexp  in  Table  6  using  Eqs.  (13)  and  (14).  The  result  derived 
from  the  values  obtained  by  Chagnes  et  al.  [31,32]  are  also  shown 
for  the  LiX  solutions  in  yBL,  in  which  <rmax  and  C*  were  calculated  in 
the  same  manner  by  the  authors  of  the  present  article  from  their 
experimental  data.  In  yBL  atx  <  1.5,  the  normalized  conductivity  y 
well  follows  the  law  of  corresponding  state,  Eq.  (18),  for  either 
solute.  In  PC,  on  the  other  hand,  the  deviation  in  the  experimental 
data  from  Eq.  (18)  is  appreciable  even  at  x  ~  1.0.  For  either  PC  or 
yBL  at  x  >  1.5,  the  normalized  conductivity  y  is  significantly 
underestimated  from  what  is  expected  by  the  cubic  root  law,  sug¬ 
gesting  its  applicable  limit.  As  indicated  by  the  ni/n2  scale  in  Fig.  3, 
the  solvent  to  solute  ratio  becomes  close  to  the  number  of  the 
solvent  molecules  that  solvate  the  solute  molecule  at  these  high 
concentrations.  For  example,  the  solvation  number  of  Li+  in  LiClC>4 
dissolved  in  ethylene  carbonate  is  determined  to  be  around  four  at 
C/mol  dm  3  <  1  [45].  The  mechanism  that  governs  the  transport 
phenomena  seems  to  be  modeled  from  the  other  viewpoint. 

The  present  study  provides  only  macroscopic  data  and  their 
empirical  treatment  which  must  be  complemented  by  microscopic 
studies  such  as  Raman  spectroscopy  to  elicit  the  molecular  images 
of  the  Na+  solvation  as  in  the  case  of  Li+.  We  also  hope  that  some 
theory  can  quantitatively  bridge  the  gap  between  macroscopic 
properties  and  microscopic  information  for  highly  concentrated 
electrolyte  solutions. 


5.  Conclusion 

While  the  partial  volume  (Vj)  of  the  NaC104  solution  is  larger 
than  that  of  the  LiClC>4  solution,  the  fluidity  (I/17)  and  the  conduc¬ 
tivity  (<r  or  A)  of  the  former  are  significantly  higher  than  those  of 
the  latter  throughout  the  concentration  up  to  1.5—2  mol  dnr3. 
Sodium-based  solutions  can  serve  as  a  good  electrolyte  for  the 
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electrochemical  devices  as  lithium-based  ones  at  least  with  respect 
to  the  conductivity.  The  maxima  in  the  specific  conductivity  (omax) 
roughly  coincide  with  those  expected  from  the  empirical  cubic  root 
law,  ffmax/C*  =  Ao/4.  The  slope  of  the  molar  conductivity  A  vs.  the 
cubic  root  of  the  molarity  C1^3  is  semi-quantitatively  interpretable 
in  the  framework  of  the  pseudolattice  model.  Deviations  from  the 
cubic  root  law  are  more  conspicuous  in  PC  than  in  yBL. 
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